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A recently proposed approximate large-n ground state solution of the Ul(ji) x Ur(u) symmetric 
linear sigma model is investigated at finite temperature. We study the coupled evaporation of two 
condensates corresponding to the symmetry breaking pattern Ul{3) x Ur(3) —> LV(2), realized by 
the ground state in certain parts of the coupling space. The region of the fluctuation induced first 
order transitions and its second order boundary is mapped out. The existence of a tricritical point 
is conjectured. 
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I. INTRODUCTION 



^ I The Lagrangian of QCD displays an approximate chiral Ul(3) x Ur(3) symmetry. Although the axial anomaly 
i reduces the approximate symmetry to SUa{3) x SU\r(3) x Uv(1), at finite temperature a fast restoration of the 
Ua(1) factor might give phenomenological relevance to investigations of the £/z(3) x £/r(3) symmetric linear sigma 
^ (— i model. The interest for investigating the symmetry breaking phase structure of the Ur,(n) x Un(n) symmetric model 
i— i at large-n, and the symmetry restoration pattern at finite temperature stems from the fact that according to the 
1 leading order renormalization group analysis, all models with n > 2 go through a first order symmetry restoration of 
^2 ■ similar characterization l], ■ The fluctuation induced first order nature [!, |j] of the transition was confirmed with 
■ nonperturbative studies Etjp, while purely perturbative methods [Tol.[Tl| or large-n approximations based on a 0(2n 2 ) 
I | symmetric saddle point 12] could enforce a region of first order transitions only by introducing the explicitly cubic 
i— ( ■ ?7a(1) anomaly term of 't Hooft. 

| It is puzzling that even the most recent lattice simulations failed to establish the region of the explicit symmetry 
* . breaking fields where the QCD phase transition displays a first order discontinuity [13l 4l5| . Predictions of the effective 
models concerning the size of this region fall in the ballpark of the lattice estimations jl6l. In particular, meson 
£N| . models solved in the limit of large flavor numbers represent a promising approach for strongly coupled effective models, 
t-H ■ when their solution is continued to n = 2, 3 [l2l. [l8j. 

A fundamental obstacle in exploring the nature of the transition with the large-n technique was the absence of any 
result for the leading order ground state, qualitatively different from the large-n behavior of the 0(2n 2 ) symmetric 
linear sigma mod el Il9l. Recently we have constructed an approximate large-n ground state for the E/t(n) X Un{n) 
^ . symmetric model |20l |21| based on the extra assumption that the heavy scalar fields can be treated at low temperatures 
as static. Therefore, in the computation of the effective potential, only pseudoscalars were included. This represents 
a consistent approximation up to temperatures where the fluctuations of the heavy scalars are excited. This feature 
is increasingly difficult to maintain near the ho = axis, where the mass gap between the scalars and pseudoscalars 
tends to zero at criticality. For this reason the findings of the present report concerning this region should be checked 
with more complete studies in the future. 

The proposed solution goes beyond the usual 1-loop effective potential, because it is built on a nontrivial saddle point 
searched in the space of a rather general set of quadratic auxiliary fields: a Lorentz-scalar Uy(n) singlet (x), a scalar 
Uv(n) vector (j/"), and a pseudoscalar Uy(n) vector (y%), introduced in [21| via Hubbard-Stratonovich transformations. 
Including the propagation of these composite fields, one effectively takes into account diagrams with more than 1 loop 
in terms of the original field variables, even when considering the 1-loop effective potential in terms of the extended 
formulation (20| . The symmetry breaking pattern ^(n) x Ur(ti) —> Uv{n— 1) investigated below is parametrized by two 
diagonal elements of the U(n) Lie algebra: v=y/ 2n 2 (vo 1 + VsT^ iag ) , where T£ iag — diag(l, 1, 1, —n + l)/y/n{n — 1) is 
the "longest" element of the Cartan subalgebra. The choice of this condensate is motivated phenomenologically, since 
it can be connected the most economically with the n = 3 case. Below, we consistently use the notation T% iag — T 8 and 
a similar notation for the corresponding U (n)-vector components. The condensate induces nonzero saddle point values 
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for the composite fields x,y®,yi. In the absence of any explicit symmetry breaking conjugate to vg = Xs/\/n, we have 
found in [21| that the configuration vo ^ 0, xs — represents the absolute minimum of the effective potential for all 
values of the quartic couplings g± . go and the external field ho conjugate to Vq, in agreement with common expectations. 

The new observation made in [2l| was the existence of a metastable local minimum even in the chiral limit, charac- 
terized by xs 7^ 0, wq 7^ 0. The corresponding vq was the same as in the ground state with the Uv(n) symmetry. The 
present investigation focuses on a part of the parameter space, where the explicit breaking hg,T s drives the newfound 
local minimum to be the true ground state at low temperatures. We determine the region of the first order transitions 
and its second order boundary. Characteristic features will be illustrated in a two-dimensional slice, where h$ and the 
0(2n 2 )-symmetric quartic coupling are fixed. Qualitative exploration is also outlined varying these two parameters. 
We speculate that at the edge of applicability of our approximation the observed features imply the existence of a 
tricritical point. 

The mixing of the composite fields with the elementary pseudoscalars results in a mass splitting within this sector. 
We found 2(n — 1) mixing heavier pseudoscalars [2l|, they are called below "kaons," since they coincide with them at 
n — 3. The (n — l) 2 light (Goldstonc) fields bear the short name "pion." There is an additional heavy U(n—1) singlet 
analogue of ?788- The 1-loop effective potential with 0(n) accuracy, depending on 2 condensates and 3 saddle-point 



coordinates (i.e. 
details in [21}): 



vq, x&^ x > Vii Ui) nas the following expression, when only the pseudoscalar contribution is included (see 



v = v cl + v„, 
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(n 2 - 2n) / ln(-;p 2 + M 2 ) + 2n / ln(-p 2 + M 2 K ) 



(1) 



One recognizes that the expression V c i contains beyond the leading terms (~ n 2 ) also subleading contributions (~ n), 
which reflect the scaling features of the fields X8,J/i,^8- The same distinction in V q reflects the multiplicities in the 
pseudoscalar sector. The saddle point equations which determine the values of the auxiliary fields are the following: 



-x = 2v 2 + T F (M 2 ) + -(xi+T F (M 2 K )-T F (M 2 ))=- i J^y° 1 



i^-yl = 2voXs-xl-l(T F (M 2 K )-T F (M 2 )), 



(2) 



where T F (n 2 ) is the finite part of a tadpole integral |20l. [21|: 
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The first term of the right-hand side is the renormalized vacuum piece at renormalization scale (fixed to m in this 
study), while the second one is the finite temperature part at temperature T. We do not show any counterterms 
explicitly; the complete renormalization program of the approximation was presented in [2l| . The masses M 2 and 
appearing in ([T]) and ([2]) are determined by the gap equations 



M 2 = M 



Ml 




(4) 



where 



M 2 



-m 



-m 2 + (gi + g 2 ) 



2v 2 +T F (M 2 ) + -(xl + T F (M 2 K ) 



T F (M 2 )) 



(5) 



The O(n ) part of ([5]) coincides with the gap equation of the 0(2n 2 ) symmetric model up to the missing scalar 
contribution, omitted by the heaviness assumption. Only the O ( 1 /n) contributions introduce intrinsic UlX Ur breaking 
features. (Note the change of the m 2 sign convention relative to |21(). The strength of all these terms is governed by 
<?2 alone, which appears in the masses figuring in the 0(1/ n) piece of this equation. This is made explicit also in ([8} 
below. Exploiting the saddle point equation of yf one obtains a gap equation for M|- which determines directly its 
dependence on vq and x§: 



Ml = M 2 - 2g 2 



X8(v - Xs) - j(T' 



\M 2 k )-T f {MD) 



(6) 
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II. APPLICATION OF THE LARGE- TV POTENTIAL AT FINITE N 



The vq 7^ 0, xs 7^ nontrivial solution originally was found in a strict large-n study by separately minimizing the 
0(n 2 ) and the 0(n) parts of the potential (JIJ [21j. Then the leading order equation describes the restoration of the 
0(2n 2 ) C Ul(ti) x Un(n) symmetry. In the subleading equations, one searched for minima with respect to xs m the 
background of vq determined from the leading 0(n 2 ) potential. In the present paper we shall consider the case n = 3, 
and search for stationary points of the sum V (both at T = and at nonzero temperatures). 



For the evaluation of the sum we use the 0(l/n) accurate values of vq, M , xs 



Ml. 



Since in the equations of xs and 



Mjj only the 0(n) part of the potential is involved, those equations are unchanged as compared to the equations used 
in our previous paper [2l|. Namely, in addition to ©, one has the following equation of state for xs after substituting 
yf from (J5J) into the derivative of V with respect to xs '■ 



= h 8 - M 2 xs - 2# 2 (>o - Xs) 



Xs(2v -Xs)-l(T F (M 2 K )-T F (M 2 )) 



-g 2X sT F (M 2 K ). 



The equation of state of Vq, including also the 0(1/ n) contributions, looks like 

= 7*0- M 2 v - ^-xs (ivoxs - xl ~ \ (T F (M 2 K ) - T F (M 2 ))^j 



(7) 



(8) 



For the realization of the strategy outlined above, Eqs. ([5]) and (|5J) were solved in the spirit of the large-n approxi- 
mation itcratively, by substituting into them 
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Ml 
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Ml 



1 

Vq = V 00 + -Vol- 



(9) 



In © and it is sufficient to replace everywhere M 2 (and M 2 ) by Mq and use vqq. These quantities (together 
with M 2 and Vox) are found from the order by order solution of Eqs. ([5]) and (|SJ) keeping the first two orders of their 
respective power series in 1/n. The values of vo and M 2 are obtained then from ©, where we put n — 3. These 
values are used together with xs and when one searches for the minima of the potential ([T]) as a function of the 
temperature. 



III. PHASE DIAGRAM FROM THE NUMERICAL SOLUTION 



Our interest in the numerical investigation was to see if the two condensates vo and v$ = Xs/V^ evaporate strongly 
correlated enough to expect a single unique transition in a certain part of the coupling space. We consider only the case 
when at T — the "nontrivial" xs minimum, not proportional to /ig represents the true ground state. The complete 
exploration of the four-dimensional parameter space (ho, hs, ffi, 32) is beyond the present scope. Since vqo depends only 
on <?i +52, we decided to fix this sum. The leading order solutions Mq ,i>oo were therefore tuned exclusively by Hq. 
We varied gi at T — in a large interval, where the approximate symmetry of the ground state is broken down to 
LV(2), if a large enough value is fixed for hs, and ho is near zero. With the choice of hs = 0.04m 3 and g\ + g% = 14.5, 
we found exclusively positive squared masses in the whole i>s interval of interest, avoiding any problem related to a 
complex effective potential 0, [23[ ■ We shall present the results for this case in some detail, followed by statements 
concerning the effect of the variation of hs and g\ +32- 

Two further observations set upper bounds on the variation of gi and ho- For gi > gimax — 14.5, the system becomes 
unstable, since in the stability region g\ > 0. Furthermore, when one increases ho at fixed gi, one encounters a value 
where the strength of hs is not sufficient anymore to drive the nontrivial solution of vs to be the absolute minimum at 
T = 0. 

The region of first order transitions restoring the approximate £V(3) symmetry was mapped out in the (g 2 ,ho) 
plane by observing a discontinuity of vs at some transition temperature and nonvanishing value of vo ■ The size of the 
discontinuity Ai>8 increases with g 2 at fixed ho- This tendency is the same as found with exact renormalization group 
studies in the n — 2 case [5j-|8|. In the opposite direction, one finds a critical value g 2 ,c(ho) below which the transition 
becomes a crossover. 

The edge of the first order transition region was localized by looking for the point where Avs vanishes. /io,c(ff2) was 
found at fixed g 2 (2.5 < g 2 < 8.5), by varying ho- Figure Q] illustrates the method. In another approach, we used 
that continuous transitions are signaled by the A singularity of an appropriate susceptibility. Remarkably, vo reacts 
sensitively to the restoration of the approximate t/y(3) symmetry becoming continuous. Namely, dvo/dho displays 
a well-expressed A singularity, although this quantity is not generically related to the Uy(2) t/y(3) change in the 
approximate symmetry pattern. This behavior reflects a rather strong coupling between the two order parameters. 
The first method locates /io,c(<72) with an accuracy 0(1O -5 ), while determination of the same curve as g 2tC (ho) was 
consistent with it. In Fig. [5] we display the points of the boundary of the first order transition region. The boundary 
which separates regions with three resp. two minima of the effective potential at T = is drawn by the dashed line. 
Correspondingly, it is possible to have the nontrivial Xs 7^ minimum as a stable or metastable state. In the latter 
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g 1= 7.5, g 2 =7.0, h a /m 3 =0.04 
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FIG. 1: The critical value of ho, c is found as the vanishing point of the discontinuity in ug when ho is varied. The curves from 
left to right display the discontinuities belonging to decreasing ho values. The critical value is fitted to ho,c ~ 0.00211m 3 . 



g 1+ g 2 =14.5, h a /m 3 =0.04 



14 



12 



6 







\ 








1 st ORDER 




• ale 




(stability) 


CROSSOVER 








(metastability) 




CROSSOVER 
(no multiple solutions) 







0.01 0.02 0.03 0.04 0.05 

h /m 3 

FIG. 2: Boundary lines of the first order transition region in the <?2 — fro plane. The quasihorizontal curve shows the second order 
boundary of the first order transition, while the quasivertical one delimits the region where h% = 0.04m 3 is sufficient to push the 
nontrivial local minimum into the global minimum of the effective potential. The dashed line only signals the boundary of the 
existence of the nontrivial (vs 7^ 0) local minimum at T = 0. 

case the phase evolution proceeds from the "trivial" vo 7^ 0, xs ~ ^8 minimum into the "fully trivial" minimum where 
wo, v s are both proportional to the explicit breaking fields. 
A high quality fit to the phase boundary of the scaling form 

ho, c = a{g 2 - bf, a = (8.70 ±0.25) • 1CT 4 , b = 8.49 ±0.01, y = 2.23 ±0.02 (10) 

shows that the phase boundary crosses the ho = axis at g 2 c « 8.5. The fit is mostly sensitive to the boundary 
points near the ho = axis, therefore, in its neighborhood the transition points were determined quite densely. The 
errors vary with the chosen points whereupon the fit was made. To obtain (fTU)) . we chose the region g 2 € [6.5,8.4]. 
Errors of the exponent should be estimated from the variation of the fitted parameters with respect to the chosen 
region of the points in the first place, not just by root- mean-square errors of a single fit. It is quite remarkable that 
the value of the critical exponent y is very close to the value 5/2, which reminds us the reciprocal of the mean field 
exponent characterizing the quark- mass/external field dep endence of various quantities (e.g. chemical potential, critical 
temperature, or couplings) near a tricritical point [24j, |25j . In the suspected tricritical point then one would expect that 
a second order line along the ho = axis starting from the origin will meet the one displayed in Fig. [2J which would 
indicate continuous transition of the two nontrivial minima into the fully trivial symmetric state at the tricritical point. 
Moreover, there is an additional interesting indication pointing to an increasingly correlated evaporation of vq and v% 
near ho = 0. 

Figure [3] shows the variation of the critical temperature T c ^(g2) along the phase boundary and the pseudocritical 
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FIG. 3: The variation of the transition temperature of vs along the phase boundary parametrized by t/2- The lower curve comes 
from the zero of the discontinuity Avg and the lupper from the maximum of dvoo/dho (which depends on only due to the 
(^-dependence of the phase boundary line) . 

temperature T c ${g-i). The latter characterizes the vq evaporation from the "trivial" into the "fully trivial" minimum, 
determined by the maximum of dvoo/dh evaluated at fto,c(<?2)- The two curves approach each other very close near 
the fig axis. It shows that the transition of vg from the nontrivial into the trivial (vg ~ hg) minimum precedes the 
transition point of uoo, an d the Vg jump is followed by a separate crossover in vqo accompanied by the further gradual 
decrease of vg. This corresponds to the sequential restoration LV(2) — > £/y(3) — > Ul(3) x Ur(3). 

With the increase of hg we observe an increase in g% c and the two curves on Fig. [3] further approach each other. 
This tendency continues monotonically until g\ c — 14.5 is reached, but the difference of the critical temperatures (AT) 
does not change sign. Further increasing hg pulls away the critical point from the ho axis along the 32 = 14.5 line 
with increasing AT. The analysis can be repeated with increasing g\ + gi values and along the ho axis. One finds 
diminishing minimum values for AT reached always at g\ = 0,<72 = <?2,c- The corresponding hg gets smaller as we 
increase g\ + 52- It is at <7i + 52 ~ 30.8 where the minimum of AT becomes zero (we have hg ~ 0.037m 3 at this point). 
Going higher in g\ + gi one even finds that the restoration of Vq precedes the evaporation of vg. This clearly shows 
that in this region the one-step recovery of the approximate full symmetry from approximate Uv(2) does occur. 

The last issue to be discussed is the validity of the assumed mass-hierarchy in the temperature range of the investiga- 
tion. The range of temperatures involved in the study does not exceed the mass scale m. For the squared mass-ratios 
in the interior of the first order region, we found that the temperature dependent behavior of the pion and kaon mass is 
compatible with neglecting the heavy scalars. However, in points near the boundary of the first order region and close 
to ho = the scalar-pseudoscalar mass difference parametrically proportional to «q decreases when the temperature 
is close to the critical value. Therefore just in the region of the conjectured tricritical behavior our results should be 
checked with an improved treatment including the scalar fluctuations, too. 

In summary, for fixed values of hg and g\ + gi we determined the boundary between the region of the crossover 
and the first order transition restoring the approximate Uy(3) symmetry and found evidence for a strongly coupled 
behavior of the condensates vo and vg near the ho — > (hg 7^ 0) limit. Tendencies observed when hg and g\ + gi are 
varied appropriately led us to a region of single step &V(2) — > Ul(3) x Ur(3) transition. 

Acknowledgements 

The authors are grateful to A. Jakovac and Zs. Szep for their careful reading of the manuscript and valuable 
suggestions. This work is supported by the Hungarian Research Fund under Contracts No. T068108 and K77534. G. 
F. is supported by the Japan Society for the Promotion of Science under ID No. P11795. 



[1] A.J. Paterson, Nucl. Phys. B 190, 188 (1981). 

[2] R.D. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 (1984). 

[3] SR. Coleman and E.J. Weinberg, Phys. Rev. D 7, 1888 (1973). 

[4] B.I. Halperin, T.C. Lubensky and S.K. Ma, Phys. Rev. Lett. 32, 292 (1974). 

[5] J. Berges, N. Tetradis and C. Wetterich, Phys. Rep. 363, 223 (2002). 

[6] J. Berges and C. Wetterich, Nucl. Phys. B 487, 675 (1997). 

[7] J. Berges, N. Tetradis and C. Wetterich, Phys. Lett. B 393, 387 (1997). 

[8] K. Fukushima, K. Kamikado and B. Klein, Phys. Rev. D 83, 116005 (2011). 



G 



[9] J.T. Lenaghan, D.H. Rischke and J. Schaffner-Bielich, Phys. Rev. D 62, 085008 (2000). 
[10] L.-H. Chan and R.W. Haymaker, Phys. Rev. D 7, 415 (1973). 
[11] C. Rosenzweig, J. Schechter and C.G. Trahern, Phys. Rev. D 21, 3388 (1980). 
[12] H. Meyer-Ortmanns and B.J. Schaefer, Phys. Rev. D 53, 6586 (1996). 

[13] F. Karsch, C. R. Allton, S. Ejiria, S. J. Hands, O. Kaczmarek, E. Laermann, and C. Schmidt, Nucl. Phys. Proc. Suppl. 129, 
614 (2004). 

[14] G. Endrodi, Z. Fodor, S.D. Katz, and K. K. Szabo, Proc. Sci. LATTICE2007, 182 (2007). 

[15] H.-T. Ding, A. Bazavov, P. Hegde, F. Karsch, S. Mukherjee, and P. Petreczky, in Proceedings of the XXIX International 

Symposium on Lattice Field Theory, Lattice 2011, Lake Tahoe, California, 2011 (to be published). 
[16] J.T. Lenaghan, Phys. Rev. D 63, 037901 (2001). 

[17] T. Herpay, A.Patkos, Zs. Szep and P. Szepfalusy, Phys. Rev. D 71, 125017 (2005). 
[18] G. Marko and Zs. Szep, Phys. Rev. B 82 065021 (2010). 
[19] L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320 (1974). 
[20] G. Fejos and A. Patkos, Phys. Rev. D 82, 045011 (2010). 
[21] G. Fejos and A. Patkos, Phys. Rev. D 84, 036001 (2011). 

[22] R.J. Rivers, Path Integral Methods in Quantum Field Theory (Cambridge Univ. Press, Cambridge, UK., 1987). 
[23] E.J. Weinberg and A. Wu, Phys. Rev. D 36, 2474 (1987). 
[24] Y. Hatta and T. Ikeda, Phys. Rev. D 67, 014028 (2003). 
[25] T. Herpay and Zs. Szep, Phys. Rev. D 74, 025008 (2006). 



